Phenotypic characterization of liver tissue heterogeneity through a next-generation 3D single-cell atlas.

Three-dimensional (3D) geometrical models are potent tools for quantifying complex tissue features and exploring structure-function relationships. However, these models are generally incomplete due to experimental limitations in acquiring multiple (> 4) fluorescent channels in thick tissue sections simultaneously. Indeed, predictive geometrical and functional models of the liver have been restricted to few tissue and cellular components, excluding important cellular populations such as hepatic stellate cells (HSCs) and Kupffer cells (KCs). Here, we combined deep-tissue immunostaining, multiphoton microscopy, deep-learning techniques, and 3D image processing to computationally expand the number of simultaneously reconstructed tissue structures. We then generated a spatial single-cell atlas of hepatic architecture (Hep3D), including all main tissue and cellular components at different stages of post-natal development in mice. We used Hep3D to quantitatively study 1) hepatic morphodynamics from early post-natal development to adulthood, and 2) the effect on the liver's overall structure when changing the hepatic environment after removing KCs. In addition to a complete description of bile canaliculi and sinusoidal network remodeling, our analysis uncovered unexpected spatiotemporal patterns of non-parenchymal cells and hepatocytes differing in size, number of nuclei, and DNA content. Surprisingly, we found that the specific depletion of KCs results in morphological changes in hepatocytes and HSCs. These findings reveal novel characteristics of liver heterogeneity and have important implications for both the structural organization of liver tissue and its function. Our next-gen 3D single-cell atlas is a powerful tool to understand liver tissue architecture, opening up avenues for in-depth investigations into tissue structure across both normal and pathological conditions.

Periphery of porcine hepatic lobes has the smallest length density of hepatic sinusoids and bile canaliculi: A stereological histological study with implications for liver biopsies.

BACKGROUND: Porcine liver is widely used in hepatologic research as a large animal model with many anatomical and physiological similarities with humans. However, only limited information on porcine liver spatial microstructure has been published, especially regarding the hepatic sinusoids and bile canaliculi. The aim of our study was to quantify the sinusoidal and bile canalicular network in healthy male and female porcine livers and to map the variability of these structures with heterogenous distribution to improve the evaluability of liver biopsy samples. METHODS: Livers from 12 healthy piglets (6 females and 6 neutered males) were sampled into 36 tissue samples per organ, representing six hepatic lobes and three different regions related to the hepatic vasculature (peripheral, paracaval and paraportal region). Histological sections were processed with a random orientation of the cutting plane. The endothelium and the bile canaliculi were stained using Ricinus communis agglutinin I lectin histochemistry. The length densities of hepatic sinusoids LV(sinusoids,liver), of bile canaliculi LV(bile canaliculi,liver) and volume fraction VV(sinusoids,liver) and surface density SV(sinusoids,liver) of sinusoids were estimated using stereological methods. The newly acquired morphometric data were compared with previously published data on density of porcine hepatocytes and fractions of connective tissue. RESULTS: The peripheral region had smallest LV(sinusoids,liver), smallest LV(bile canaliculi,liver) and greatest VV(sinusoids,liver). The six hepatic lobes had statistically comparable length densities of both sinusoids and bile canaliculi, but the left lateral lobe had smallest VV(sinusoids,liver). Regions with greater LV(sinusoids,liver) had also greater LV(bile canaliculi,liver) and SV(sinusoids,liver) and were accompanied by greater density of smaller hepatocytes. Regions with smaller LV(sinusoids,liver) and LV(bile canaliculi,liver) contained a greater fraction of interlobular connective tissue. CONCLUSIONS: The length density of hepatic sinusoids is smaller in the peripheral regions of the porcine liver than in other regions related to the hepatic vasculature - paracaval and paraportal regions, and smaller in castrated males than in females. Greater length density of liver sinusoids was linked with greater local density of bile canaliculi, with local increase in the density of smaller hepatocytes and, simultaneously, with smaller fractions of hepatic connective tissue. The intrahepatic and inter-sexual variability of the porcine liver morphology needs to be taken into account when designing and interpreting experiments involving the histological quantification of the microvascular network. The complete primary morphometric data describing the distribution of morphometric parameters within porcine liver were made available in a form facilitating the power analysis to justify the minimal number of tissue samples or animals required when designing further histological evaluation studies. The macroscopic map of microvessels and bile canaliculi variability facilitates their assessment in liver biopsies in the pig.

Induction of open-form bile canaliculus formation by hepatocytes for evaluation of biliary drug excretion.

Biliary excretion is a major drug elimination pathway that affects their efficacy and safety. The currently available in vitro sandwich-cultured hepatocyte method is cumbersome because drugs accumulate in the closed bile canalicular lumen formed between hepatocytes and their amounts cannot be mealsured directly. This study proposes a hepatocyte culture model for the rapid evaluation of drug biliary excretion using permeation assays. When hepatocytes are cultured on a permeable support coated with the cell adhesion protein claudins, an open-form bile canalicular lumen is formed at the surface of the permeable support. Upon application to the basolateral (blood) side, drugs appear on the bile canalicular side. The biliary excretion clearance of several drugs, as estimated from the obtained permeabilities, correlates well with the reported in vivo biliary excretion clearance in humans. Thus, the established model is useful for applications in the efficient evaluation of biliary excretion during drug discovery and development.

Apical bulkheads accumulate as adaptive response to impaired bile flow in liver disease.

Hepatocytes form bile canaliculi that dynamically respond to the signalling activity of bile acids and bile flow. Little is known about their responses to intraluminal pressure. During embryonic development, hepatocytes assemble apical bulkheads that increase the canalicular resistance to intraluminal pressure. Here, we investigate whether they also protect bile canaliculi against elevated pressure upon impaired bile flow in adult liver. Apical bulkheads accumulate upon bile flow obstruction in mouse models and patients with primary sclerosing cholangitis (PSC). Their loss under these conditions leads to abnormally dilated canaliculi, resembling liver cell rosettes described in other hepatic diseases. 3D reconstruction reveals that these structures are sections of cysts and tubes formed by hepatocytes. Mathematical modelling establishes that they positively correlate with canalicular pressure and occur in early PSC stages. Using primary hepatocytes and 3D organoids, we demonstrate that excessive canalicular pressure causes the loss of apical bulkheads and formation of rosettes. Our results suggest that apical bulkheads are a protective mechanism of hepatocytes against impaired bile flow, highlighting the role of canalicular pressure in liver diseases.

A modular analysis of bile canalicular function and its implications for cholestasis.

Hepatocytes produce bile components and secrete them into a lumen, known as a bile canaliculus, that is formed by the apical membranes of adjoining hepatocytes. Bile canaliculi merge to form tubular structures that subsequently connect to the canal of Hering and larger intra- and extrahepatic bile ducts formed by cholangiocytes, which modify bile and enable flow through the small intestine. The major functional requirements for bile canaliculi are the maintenance of canalicular shape to preserve the blood-bile barrier and regulation of bile flow. These functional requirements are mediated by functional modules, primarily transporters, the cytoskeleton, cell-cell junctions, and mechanosensing proteins. I propose here that bile canaliculi behave as robust machines whereby the functional modules act in a coordinated manner to perform the multistep task of maintaining canalicular shape and bile flow. Cholestasis, the general term for aberrant bile flow, stems from drug/toxin-induced or genetic dysregulation of one or more of the protein components in the functional modules. Here, I discuss the interactions between components of the various functional modules in bile canaliculi and describe how these functional modules regulate canalicular morphology and function. I use this framework to provide a perspective on recent studies of bile canalicular dynamics.

Construction of extended and functional bile canaliculi using long-term sandwich-cultured cryopreserved human hepatocytes and the application of hepatocytes for predicting the biliary excretion of pharmaceutical and food-related compounds.

The biliary excretion of pharmaceutical and food-related compounds is an important factor for assessing pharmacokinetics and toxicities in humans, and a highly predictive in vitro method for human biliary excretion is required. We have developed a simple in vitro culture method for generating extended and functional bile canaliculi using cryopreserved human hepatocytes. We evaluated the uptake of compounds by hepatocytes and bile canaliculi, and the biliary excretion index (BEI) was calculated. After 21 days of culture, the presence of extended and functional bile canaliculi was confirmed by the uptake of two fluorescent substrates. Positive BEIs were observed for taurocholic acid-d4, rosuvastatin, pitavastatin, pravastatin, valsartan, olmesartan, and topotecan (reported biliary-excreted compounds in humans), but no difference in BEI was observed for salicylic acid (a nonbiliary-excreted compound). Furthermore, 8 of 21 food-related compounds with specific structures and reported biliary transporter involvement exhibited positive BEIs. The developed in vitro system was characterized by functional bile canaliculus-like structures, and it could be applied to the prediction of the biliary excretion of pharmaceutical and food-related compounds.

Hepatocyte apical bulkheads provide a mechanical means to oppose bile pressure.

Hepatocytes grow their apical surfaces anisotropically to generate a 3D network of bile canaliculi (BC). BC elongation is ensured by apical bulkheads, membrane extensions that traverse the lumen and connect juxtaposed hepatocytes. We hypothesize that apical bulkheads are mechanical elements that shape the BC lumen in liver development but also counteract elevated biliary pressure. Here, by resolving their structure using STED microscopy, we found that they are sealed by tight junction loops, connected by adherens junctions, and contain contractile actomyosin, characteristics of mechanical function. Apical bulkheads persist at high pressure upon microinjection of fluid into the BC lumen, and laser ablation demonstrated that they are under tension. A mechanical model based on ablation results revealed that apical bulkheads double the pressure BC can hold. Apical bulkhead frequency anticorrelates with BC connectivity during mouse liver development, consistent with predicted changes in biliary pressure. Our findings demonstrate that apical bulkheads are load-bearing mechanical elements that could protect the BC network against elevated pressure.

The morphology, functionality, and longevity of a novel all human hepatic cell-based tri-culture system.

There remains a significant need for a convenient, phenotypically stable long-term culture platform for primary human hepatocytes (PHHs) for use in pharmacological and toxicological applications. Conventional in vitro models are often inconvenient, burdensome to use, and unable to support a multitude of donor lots or maintain PHH structural and functional properties over extended time. To address these limitations, an all-human cell-based hepatic tri-culture system (HTCS) has been developed comprised of frozen vials of PHHs and feeder cells. Qualified PHHs exhibited healthy morphological characteristics for ≥30 days. Extensive anastomosing networks of bile canaliculi with tight and gap junctions were established early and remained stable and functional throughout the culture period. After 5 culture days, albumin, urea, and basal Phase 1 and Phase 2 metabolic functions were stable for at least 2 weeks and significantly higher in the HTCS PHHs compared to sandwich monoculture PHHs. Induction of CYP functional activity by prototypical receptor agonists was stable after 4 days for at least 2 weeks. Gene expression of Alb and various CYPs in the HTCS PHHs was significantly higher compared to sandwich monoculture PHHs. The HTCS represents a convenient, phenotypically stable, all-human PHH culture platform for pharmacological and toxicological applications.

Induction of Bile Canaliculi-Forming Hepatocytes from Human Pluripotent Stem Cells.

Cell polarity and formation of bile canaliculi can be achieved in hepatocytes which are generated from patient-derived induced pluripotent stem cells. This allows for the study of endogenous mutant proteins, patient-specific pathogenesis, and drug responses for diseases where hepatocyte polarity and bile canaliculi play a key role. Here, we describe a step-by-step protocol for the generation of bile canaliculi-forming hepatocytes from induced pluripotent stem cells and their evaluation.

Construction of a culture protocol for functional bile canaliculi formation to apply human iPS cell-derived hepatocytes for cholestasis evaluation.

Cholestatic toxicity causes the failure of pharmaceutical agents during drug development and, thus, should be identified at an early stage of drug discovery and development. The formation of functional bile canaliculi in human hepatocytes is required for in vitro cholestasis toxicity tests conducted during the early stage of drug development. In this study, we investigated the culture conditions required for the formation of bile canaliculi using human-induced pluripotent stem cell-derived hepatocytes (hiPSC-Heps). When hiPSC-Heps were sandwich-cultured under the condition we established, extended bile canaliculi were formed on the whole well surfaces. Biliary efflux transporters were localized in the formed bile canaliculi structures which had junctional complexes. After the model substrates of the biliary efflux transporters were taken up into cells, their subsequent excretion into the bile canaliculi was observed and was found to be impeded by each inhibitor of the biliary efflux transporter. These findings suggest that bile canaliculi have transporter-specific bile excretion abilities. We will continue to study the application of this culture protocol to cell-based cholestasis assay system. As a result, the culture protocol could lead to a highly predictable, robust cell-based cholestasis assay system because it forms functional bile canaliculi reproducibly and efficiently.

Quantitative understanding of HepaRG cells during drug-induced intrahepatic cholestasis through changes in bile canaliculi dynamics.

An understanding of the quantitative relationship between bile canaliculus (BC) dynamics and the disruption of tight junctions (TJs) during drug-induced intrahepatic cholestasis may lead to new strategies aimed at drug development and toxicity testing. To investigate the relationship between BC dynamics and TJ disruption, we retrospectively analyzed the extent of TJ disruption in response to changes in the dynamics of BCs cultured with entacapone (ENT). Three hours after adding ENT, the ZO-1-negative BC surface area ratio became significantly higher (4.1-fold) than those of ZO-1-positive BCs. Based on these data, we calculated slopes of surface area changes, m, of each ZO-1-positive and ZO-1-negative BC. BCs with m ≤ 15 that fell within the 95% confidence interval of ZO-1-positive BCs were defined as ZO-1-positive. To validate this method, we compared the frequency of ZO-1-positive BCs, FZ , with that of BCs with m ≤ 15, FT , in culture using drugs that regulate TJ, or induce intrahepatic cholestasis. FT values were correlated with FZ under all culture conditions (R2  = .99). Our results indicate that the magnitude of BC surface area changes is a factor affecting TJ disruption, suggesting that maintaining TJ integrity by slowing BC dilation inhibits cell death.

HepG2-NIAS cells, a new subline of HepG2 cells that can enhance not only CYP3A4 activity but also expression of drug transporters and form bile canaliculus-like networks by the oxygenation culture via a collagen vitrigel membrane.

We reported the enhanced liver-specific function and structure of HepG2 cells by the oxygenation culture via a collagen vitrigel membrane (CVM). The cells were conditioned in our laboratory for a long period, so their characteristics may change from the original HepG2 cells registered in RIKEN cell bank (RCB) with the number of 1648 (HepG2-RCB1648 cells). We named the conditioned HepG2-RCB1648 cells in our laboratory as HepG2-NIAS cells. Here, we clarified the features of HepG2 cells with three different culture histories by analyzing their morphology and viability, CYP3A4 activity, the potential to form bile canaliculus-like structures, and the expression of drug transporters. On plastic, HepG2-NIAS cells grew as a monolayer without the formation of large aggregates involving dead cells that were observed in HepG2-RCB1648 cells and HepG2-RCB1886 cells. In the oxygenation culture via a CVM, the CYP3A4 activity of HepG2-NIAS cells increased to almost half level in direct comparison to that of differentiated HepaRG cells cultured on a collagen-coated plate; however, that of HepG2-RCB1648 cells and HepG2-RCB1886 cells was almost not detected. HepG2-NIAS cells formed bile canaliculus-like networks in which fluorescein was accumulated after the exposure of fluorescein diacetate, although HepG2-RCB1648 cells and HepG2-RCB1886 cells did not possess the potential. Also, immunohistological observations revealed that HepG2-NIAS cells remarkably enhanced the expression of drug transporters, NTCP, OATP1B1, OATP1B3, BSEP, MDR1, MRP2, and BCRP. These results suggest that HepG2-NIAS cells are a new subline of HepG2 cells useful for drug development studies. HepG2-NIAS cells were registered in RCB with the number of 4679.

New Pharmacokinetic Parameters of Imaging Substrates Quantified from Rat Liver Compartments.

Hepatobiliary imaging is increasingly used by pharmacologists to quantify liver concentrations of transporter-dependent drugs. However, liver imaging does not quantify concentrations in extracellular space, hepatocytes, and bile canaliculi. Our study compared the compartmental distribution of two hepatobiliary substrates gadobenate dimeglumine [BOPTA; 0.08 liver extraction ratio (ER)] and mebrofenin (MEB; 0.93 ER) in a model of perfused rat liver. A gamma counter placed over livers measured liver concentrations. Livers were preperfused with gadopentetate dimeglumine to measure extracellular concentrations. Concentrations coming from bile canaliculi and hepatocytes were calculated. Transporter activities were assessed by concentration ratios between compartments and pharmacokinetic parameters that describe the accumulation and decay profiles of hepatocyte concentrations. The high liver concentrations of MEB relied mainly on hepatocyte and bile canaliculi concentrations. In contrast, the three compartments contributed to the low liver concentrations obtained during BOPTA perfusion. Nonlinear regression analysis of substrate accumulation in hepatocytes revealed that cellular efflux is measurable ∼4 minutes after the start of perfusion. The hepatocyte-to-extracellular concentration ratio measured at this time point was much higher during MEB perfusion. BOPTA transport by multidrug resistance associated protein 2 induced an aquaporin-mediated water transport, whereas MEB transport did not. BOPTA clearance from hepatocytes to bile canaliculi was higher than MEB clearance. MEB did not efflux back to sinusoids, whereas BOPTA basolateral efflux contributed to the decrease in hepatocyte concentrations. In conclusion, our ex vivo model quantifies substrate compartmental distribution and transport across hepatocyte membranes and provides an additional understanding of substrate distribution in the liver. SIGNIFICANCE STATEMENT: When transporter-dependent drugs target hepatocytes, cellular concentrations are important to investigate. Low concentrations on cellular targets impair drug therapeutic effects, whereas excessive hepatocyte concentrations may induce cellular toxicity. With a gamma counter placed over rat perfused livers, we measured substrate concentrations in the extracellular space, hepatocytes, and bile canaliculi. Transport across hepatocyte membranes was calculated. The study provides an additional understanding of substrate distribution in the liver.

Steatosis Alters the Activity of Hepatocyte Membrane Transporters in Obese Rats.

Fat accumulation (steatosis) in ballooned hepatocytes alters the expression of membrane transporters in Zucker fatty (fa/fa) rats. The aim of the study was to quantify the functions of these transporters and their impact on hepatocyte concentrations using a clinical hepatobiliary contrast agent (Gadobenate dimeglumine, BOPTA) for liver imaging. In isolated and perfused rat livers, we quantified BOPTA accumulation and decay profiles in fa/+ (normal) and fa/fa hepatocytes by placing a gamma counter over livers. Profiles of BOPTA accumulation and decay in hepatocytes were analysed with nonlinear regressions to characterise BOPTA influx and efflux across hepatocyte transporters. At the end of the accumulation period, BOPTA hepatocyte concentrations and influx clearances were not significantly different in fa/+ and fa/fa livers. In contrast, bile clearance was significantly lower in fatty hepatocytes while efflux clearance back to sinusoids compensated the low efflux into canaliculi. The time when BOPTA cellular efflux impacts the accumulation profile of hepatocyte concentrations was slightly delayed (2 min) by steatosis, anticipating a delayed emptying of hepatocytes. The experimental model is useful for quantifying the functions of hepatocyte transporters in liver diseases.

Rab35 controls formation of luminal projections required for bile canalicular morphogenesis.

Hepatocytes display a unique biaxial polarity with shared apical luminal connections between adjacent hepatocytes that merge into a network of bile canaliculi. Belicova et al. (2021. J. Cell Biol.https://doi.org/10.1083/jcb.202103003) discovered that hepatocyte apical membranes generate Rab35-dependent extensions that traverse the lumen and are essential for bile canalicular formation and maintenance.

Anisotropic expansion of hepatocyte lumina enforced by apical bulkheads.

Lumen morphogenesis results from the interplay between molecular pathways and mechanical forces. In several organs, epithelial cells share their apical surfaces to form a tubular lumen. In the liver, however, hepatocytes share the apical surface only between adjacent cells and form narrow lumina that grow anisotropically, generating a 3D network of bile canaliculi (BC). Here, by studying lumenogenesis in differentiating mouse hepatoblasts in vitro, we discovered that adjacent hepatocytes assemble a pattern of specific extensions of the apical membrane traversing the lumen and ensuring its anisotropic expansion. These previously unrecognized structures form a pattern, reminiscent of the bulkheads of boats, also present in the developing and adult liver. Silencing of Rab35 resulted in loss of apical bulkheads and lumen anisotropy, leading to cyst formation. Strikingly, we could reengineer hepatocyte polarity in embryonic liver tissue, converting BC into epithelial tubes. Our results suggest that apical bulkheads are cell-intrinsic anisotropic mechanical elements that determine the elongation of BC during liver tissue morphogenesis.

Molecular Regulation of Canalicular ABC Transporters.

The ATP-binding cassette (ABC) transporters expressed at the canalicular membrane of hepatocytes mediate the secretion of several compounds into the bile canaliculi and therefore play a key role in bile secretion. Among these transporters, ABCB11 secretes bile acids, ABCB4 translocates phosphatidylcholine and ABCG5/G8 is responsible for cholesterol secretion, while ABCB1 and ABCC2 transport a variety of drugs and other compounds. The dysfunction of these transporters leads to severe, rare, evolutionary biliary diseases. The development of new therapies for patients with these diseases requires a deep understanding of the biology of these transporters. In this review, we report the current knowledge regarding the regulation of canalicular ABC transporters' folding, trafficking, membrane stability and function, and we highlight the role of molecular partners in these regulating mechanisms.

Bile formation and secretion: An update.

Bile formation is a fundamental physiological process that is vital to the survival of all vertebrates. However, little was known about the mechanisms of this secretion until after World War II. Initial studies involved classic physiologic studies in animal models and humans, which progressed to include studies in isolated cells and membrane vesicles. The advent of molecular biology then led to the identification of specific transport systems that are the determinants of this secretion. Progress in this field was reviewed in the American Physiologic Society's series on "Comprehensive Physiology" in 2013. Herein, we provide an in-depth update of progress since that time.

Protective Functions of ZO-2/Tjp2 Expressed in Hepatocytes and Cholangiocytes Against Liver Injury and Cholestasis.

BACKGROUND & AIMS: Liver tight junctions (TJs) establish tissue barriers that isolate bile from the blood circulation. TJP2/ZO-2-inactivating mutations cause progressive cholestatic liver disease in humans. Because the underlying mechanisms remain elusive, we characterized mice with liver-specific inactivation of Tjp2. METHODS: Tjp2 was deleted in hepatocytes, cholangiocytes, or both. Effects on the liver were assessed by biochemical analyses of plasma, liver, and bile and by electron microscopy, histology, and immunostaining. TJ barrier permeability was evaluated using fluorescein isothiocyanate-dextran (4 kDa). Cholic acid (CA) diet was used to assess susceptibility to liver injury. RESULTS: Liver-specific deletion of Tjp2 resulted in lower Cldn1 protein levels, minor changes to the TJ, dilated canaliculi, lower microvilli density, and aberrant radixin and bile salt export pump (BSEP) distribution, without an overt increase in TJ permeability. Hepatic Tjp2-defcient mice presented with mild progressive cholestasis with lower expression levels of bile acid transporter Abcb11/Bsep and detoxification enzyme Cyp2b10. A CA diet tolerated by control mice caused severe cholestasis and liver necrosis in Tjp2-deficient animals. 1,4-Bis[2-(3,5-dichloropyridyloxy)]benzene ameliorated CA-induced injury by enhancing Cyp2b10 expression, and ursodeoxycholic acid provided partial improvement. Inactivating Tjp2 separately in hepatocytes or cholangiocytes showed only mild CA-induced liver injury. CONCLUSION: Tjp2 is required for normal cortical distribution of radixin, canalicular volume regulation, and microvilli density. Its inactivation deregulated expression of Cldn1 and key bile acid transporters and detoxification enzymes. The mice provide a novel animal model for cholestatic liver disease caused by TJP2-inactivating mutations in humans.